Orbital Dzyaloshinskii-Moriya Exchange Interaction 
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Superexchange calculation is performed for multi-orbital band models with broken inversion sym- 
metry. Orbital-changing hopping terms allowed by the symmetry breaking electric field lead to a 
new kind of orbital exchange term closely resembling the Dzyaloshinskii-Moriya spin exchange. The 
final superexchange Hamiltonian in two dimensions is expressed in terms of Gell-Mann matrices 
that break down as three spin operators and five nematic (quadrupole) order operators. Mean-field 
phase diagram exhibits a rich structure including anti-ferro-orbital, ferro-orbital, and both single 
and multiple spiral-orbital phases. Superexchange calculation for spinful, two-orbital model also 
strongly suggests a robust chiral-orbital order in the ground state. 
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Strong on-site repulsion transforms the Bloch bands 
into an insulator where the residual low-energy dy- 
namics are the spin degrees of freedom interacting 
with each other via the superexchange mechanism [1]. 
For spin-orbit-coupled bands, the spin-flip hopping pro- 
cesses result in a new type of spin exchange called the 
Dzyaloshinskii-Moriya (DM) interaction [2] under the 
superexchange process. Symmetry-wise, local inversion 
symmetry breaking such as the bond distortion for a pair 
of adjacent magnetic orbitals, in addition to the spin- 
orbit interaction (SOI), is the pre-requisite for the DM 
interaction to make its appearance in a given system. 
The ordered magnetic ground state is modified from the 
simple coUinear structure as a result of the DM exchange. 

Inversion symmetry breaking (ISB) on the global scale 
takes place for surfaces and interfaces and affects the 
band structure directly with new effects such as the 
Rashba interaction [3]. The situation was recently re- 
viewed carefully in Refs. [4, 5] where it was shown that 
the symmetry-breaking electric field along the z-direction 
modifies the band structure within the a:?/-plane by al- 
lowing previously forbidden hopping processes. Exam- 
ples are Px{y) ^ Pz orbital hopping in the p-band, and 
dxy dzx orbital hopping in the i2g-band. It was fur- 
ther shown [5] that the new hopping terms arising from 
ISB can be cast in the form — TX^k'^k-'-' ' ^ •S)^k 
around the T (k = 0) point, where ^k is the collection 
of, say, p-orbital operators (Xk, Vk, -^k)"^ in momentum 
coordinates k, 7 is a parameter measuring the degree 
of ISB, and L is the spin-1 orbital angular momentum 
(0AM) operator. Just from the form of the new interac- 
tion it is clear that each band will carry polarized 0AM 
proportional to /i(k x z) with the respective helicities 
h = +1,0, —1 [5]. The enlargement of effective spin size 
from 1/2 (as in electrons' spin) to 1 (as in degenerate p- 
orbital bands) results in the appearance of the third band 
that remains unpolarized. The chiral structure of the 
0AM, dubbed the "orbital Rashba effect", occurs even 



in the complete absence of SOI and has been confirmed 
by circular dichroism ARPES in the weak-SOI material, 
Cu [6]. 

The new hopping processes allowed by ISB are the or- 
bital analogues of spin-flip hoppings in spin-orbit-coupled 
bands. Therefore, the two necessary conditions for the 
emergence of spin-DM interaction - ISB and SOI - are au- 
tomatically fulfilled when a symmetry-breaking electric 
field acts perpendicular to the two-dimensional surface. 
The purpose of this paper is to review this situation care- 
fully in the limit of strong on-site interaction regime to 
ask if an orbital analogue of spin-DM interaction exists. 

We address this question in two complementary ways. 
First, the spin-orbit effect will be ignored in favor of a 
spinless multi-orbital band model with strong on-site re- 
pulsion. The proposed model may find useful applica- 
tions in spinless fermion lattice or in the bosonic cold- 
atom optical lattice problem as several recent papers on 
this subject can testify [7, 8]. The orbital version of DM 
interaction indeed arises from superexchange-type cal- 
culation, but with complications from the three-orbital 
rather than the two-spin degrees of freedom present at 
each site. The effective spin-1 model involves not only 
the usual spin-spin exchange among the spin-1 operators, 
but also the interaction among the nematic operators. 
As noted recently in the cold-atom context [8], the whole 
low-energy Hamiltonian is expressed in terms of the eight 
Gell-Mann matrices, which we further demonstrate to be 
decomposed rather nicely into a set of three spin matrices 
and the other five representing the quadrupole (nematic) 
order. 

Assuming three p-orbital states (of either hard-core 
bosons or spinless fermions) at each site, we write down 
a square lattice Hamiltonian with nearest-neighbor hop- 
ping Ht = J2iiHiA+x + Hi^i+y), 
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+^{XjZi+s - ZjXi+si) + h.c. 

Hi,i + y = taY^YiJ^y + t}j[XlXiJ^y + Z^Zi^y] 

+ ^{Y^Z, + y-ZlY, + y)+h.C. . (1) 

Two hopping integrals ta and ti, are given for a- and tt- 
bonding orbitals, respectively. Inter-orbital hopping be- 
comes possible when the ISB parameter 7 is nonzero [4, 
5]. All hopping parameters are real due to the assumed 
time-reversal invariance. Three-component spinor can be 
formed, = {Xi Yi Z^) , representing the p^-, Py-, 
and p^-orbitals at the site i. To reduce the vast com- 



plexity of the supercxchange calculation to a manage- 
able minimum we adopt a simplified multi-orbital Hub- 
bard interaction Hu = UJ2i''^i{''^i ~ 1)) where rij = 
XjXi +YfYi + zjZi is the total particle number at the 
site i. 

Superexchange calculation at 1/3-filling (one parti- 
cle per each site) can proceed with the Hamiltonian 
H = Ht + Hjj- Without much loss of generality and 
to make the final outcome of the calculation tractable 
to follow (even to write down!), we had to further drop 
the weaker tt- hopping term t^. The exchange Hamilto- 
nian thus obtained, writing U = 1 and ta = 2i, reads 
'H = Hx + • 



i i 

+ 7^ X] [''^i ~ ^i^i+x - ^i^+x + ii^i? - ^){[K+x? - 1) > 

i 

Hy = 2^ ^(Af + [A^]^ - + [A^+,]^ - 1) + 2i7X: [(^i + i^i]")^Uy " >^i>^+y + I^i+s]') 

i i 

+ 7^ E [^1^1+5 - >^>^+y - >^>^+y + ([^r]' - m^iy? - 1) 

i 

I 



(2) 



Here A" arc the Gell-Mann matrices, a = 1 • • • 8, and 
A"^ = diag(l, 0,-1), A^ = diag(0, 1, —1) are combinations 
of A^, A* and the unit matrix. 

Although the Hamiltonian obtained above appears 
very complex, Hx or Hy examined in isolation has a sim- 
ple interpretation. Note that Hx only involves the three 
operators (A^, A^, A^) that, in common, have zero matrix 
elements for the pj^-orbital state and reduce to the three 
Pauli matrices r = (r"^, r^, r^) within the {px,Pz)-orhital 
subspace. In particular the part of Hx involving the first 
power of 7 can be re- written as 



2^7 Et' 



T- T' 



i-\-x 



i-\-x. 



= 2t7^t/- [Ti X Ti+s], (3) 



2^7 y ■ [fii X Hi+y]. Such Pauli matrix description (ef- 
fective spin- 1/2 model) must give way to the full Gell- 
Mann matrix formalism shown in Eq. (2) once Hx and 
Hy are combined in the two-dimensional lattice. A sim- 
ilar Gell-Mann matrix expression appeared in the low- 
energy theory of hard-core three-component bosons con- 
fined in the optical lattice [8] where, however, the influ- 
ence of ISB on the superexchange process was not exam- 
ined. 

We can further organize the Gell-Mann matrix-based 
Hamiltonian of Eq. (2) with the following observation. 
First we define, formally, the 5=1 spin operator S = 
(5^, 5*2', 5^) = (A^-A^A2) satisfying (5°)^-,. = -iSap^. 
Then the remaining five Gell-Mann matrices are gener- 
ated as 



which is exactly the orbital analogue of the DM spin ex- 
change, or "orbital DM" (ODM) exchange. The real- 
valued transition amplitudes obtained from the superex- 
change process excludes the imaginary Ty operator, per- 
mitting y • [tj X Tj+j] as the only permissible DM in- 
teraction. The 7^-term reduces to the usual anisotropic 
spin exchange ~ [TfTj — t^Tj — t^t^] in the two-orbital 
limit. Analogously, H^ reduces to the two-component 
(pj,.p^)-orbital model with (A^,A^,A'') acting as an- 
other set of Pauli matrices /it with the DM interaction. 





= s'^sy + sys"" = 






= sys' + s'-sy = - 
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= [5"]2 - = - 


-A^ 






'-[5^f) = A^ (4) 



which are precisely the five nematic operators often em- 
ployed in the discussion of quantum nematic order in 
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spin-1 models [9, 10]. Our Hamiltonian (2) thus embod- 
ies complex exchange interactions among (pseudo-)spin 
and nematic order parameters. 

We perform the mean-field analysis of the phase dia- 
gram based on Eq. (2), by making replacement of the 
operator Af with its average, (Af), in the Hamiltonian. 
Assuming site-factorized wave function 



(a) 



(5) 



with the three complex coefficients satisfying |AfP -I- 
l^fP + l^iP ~ 1' many-body wave function be- 
comes the direct product — In this scheme 
one can choose real without loss of generality and 
parameterize the coefficients generally as Af = coscti, 



A^ 



sin Qfi cos /3i , Al = e ' sina^sin/Si. The mean- 



field Hamiltonian being a function of the four angles 
(ai, l3i,"fi,Si) per site can be minimized by the Monte 
Carlo (MC) annealing method. We enrich the phase dia- 
gram of the model (2) by introducing one more parameter 
A to replace t-f tAj. Varying the strength of A in- 
dependently of 7 permits more control of the strength of 
orbital DM exchange over the other interactions in the 
model. The zero-temperature phase diagram spanning 
(7, A) is shown in Fig. 1. Everywhere we find (S^) = 0. 
Recall (iS")^^ = —iSap-y is pure imaginary, while our 
numerical MC search only turned up real-valued wave 
functions, "/i = Si = 0. The remaining six Gell-Mann ma- 
trices, including the unit matrix, can be organized into 
two groups Tf = (A4,Af,Af) and Tf = (Af,A[,Af) as 
was previously considered for Ji^ and "Hy, respectively. 
The ground state mean-field configuration was Fourier 



analyzed TV 



X, y) in search of 



modulated structures with long periods that often ap- 
pears with the DM interaction. 

When 7 = the Hamiltonian Hx + Hy is expressed 
in terms of two commuting matrices, (Af,Aj'), with the 
ground state given by alternate occupations of Px and 
Py orbitals on the square lattice. This phase, called the 
antiferro-orbital (AFO) state, dominates the small-7 re- 
gion of the phase diagram. For A small and 7 increasing 
beyond a critical value, one finds a first-order transition 
into a ferro-orbital (FO) state with p^-orbital occupation 
at every site. Inhabiting the large-A, intcrmediate-7 re- 
gion is the (11) helical (H) phase that we found to be 
well described by 



(Tf ) = (sin2[k • r, + ao] sin[k • r,-|-/3n], 0, 

cos^[k • Yi + ao] ^ sin^[k • rj + ag] sin^[k • r^-f /3o] 
(Tf)= (sin2[k • r, + ao] sin2[k • r, + /3o], 0, 

sin^[k-rj-|-ao]cos2[k-r^-|-;3o]), (6) 
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FIG. 1: (color online) (a) Zero-temperature phase diagram 
of the orbital model in Eq. (2) with t = 1. Phase bound- 
aries are obtained on the basis of extensive MC simulation 
and variational energy calculation. Orbital configurations are 
abbreviated as AFO (antiferro-orbital), FO (ferro-orbital), H 
(helical), and MH (multiple helical). Bragg patterns for each 
phase are schematically shown with dotted hne denoting the 
first Brillouin zone boundary. In H, the period increases con- 
tinuously as 7 becomes larger. Bragg spots at (fc, k) and 
{3k, 3k) are both present in H. Real space spin configurations 
in terms of (Tf ) for (b) H and (c) MH phases are depicted. 



k — {k,k). Variational calculation of the minimum en- 
ergy with respect to k and the relative phase angle ao— /3o 
confirms that the period 27r/fc is increased with 7 in 
A > 1 region from 2 to > 6 before it is supplanted by 
another intricate phase taking place between H and FO. 
The new phase, indicated as MH (multiple helical) in Fig. 
1) is constructed as the equal- weight superposition of two 
helices with ki = ±(7r/2,7r/2) and k2 = ±(7r/2, -7r/2), 
as well as two asymmetric peaks at (tt, 0) and (0, tt). 

The problem of including the physical spin in the 
three-orbital model proved too complicated, but there 
is an alternative to study much of the same physics, 
hopefully, in a more tractable two-orbital model. Write 
down the two orbitals carrying spin a as Ui^a- and 6^.0-, 

and the spinor V'^o- = ( f*'"^ ) . The two-orbital spin 



model features hopping among the same orbitals with 
amplitudes ta and tt, respectively, and the inter-orbital 
hopping proportional to the ISB parameter 7: Ht = 
J2{ij),cr '^i,<7 (ti + ^I'^v + ^2Tz) i'j,a- The Paufi matrix r is 
introduced to span the 2x2 matrix and ti = {ta+h)/2, 
^2 = {ta — tb)/'2. Time-reversal symmetry is preserved 
for this model. Implementing separate unitary rotations 
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for ^Pi and Vj, V'i e'^^f/^Vii, ^ e-^^-'y/^j [11], if* 
becomes 

II (ti'lai'j,<r+t2i'l^T,i^j,a + h.c}j (7) 

for each adjacent pair (ij), and t(cos^,sin^) = (ti,7). 

Superexchange calculation is most conveniently carried 
out in the new basis, assuming the quarter-filling (one 



P 

Qij 

Rij 

with A = t^/{U-3JH), B = e/{U-2JH), and C = t^/U. 
The Hund's exchange Jr is included in the present cal- 
culation. All the orbital operators fi and fj arc writ- 
ten in the transformed basis ij^i and ipj^ and to recover 
the form in the original orbital basis one needs the ro- 
tation ff = rf cos 9 — rf sin 9, f f = rf cos 9 + rf sin 6*, 
and ff = rf cos 9 + r? sin 9, fj = rj cos 9 — rj sin 0. As 
a result fi • fj and rf in Eq. (9) become 

cos 26 Ti ■Tj+2 sin^ 9 t^tJ + sin 26 {y ■ n x r, ), 
cos^ 9 T^T^ - sin^ 6* rf r/ ^ sin 261 {y ■ Ti x r,), (10) 

respectively. There is again the orbital version of the 
DM interaction, y-Ti xtj, with strength proportional to 
sin 26 ~ 7. 

A simpler form follows from taking Jh = 0, 

2t 

+ (^|) (tTi.e7, + l)(2fff/-f,-f,- + l), (11) 

where the overall exchange energy J = /2U is taken 
to unity. The first term possesses the well-known SU(4) 
symmetry [12], broken down to SU(2)xU(l) by terms 
of order t^/t and {t^/tY in the above. Microscopically 
nonzero t2 arises from diff'erent hopping amplitudes of 
CT-bonding and 7r-bonding orbitals on the lattice. One 
can see the qualitative effects of the symmetry-breaking 
terms within the mean-field picture. The interaction lin- 
ear in t2/t favors a pair of adjacent spins to be anti- 
parallel, (cTj • (jj) — 1 < 0, therefore the orbital opera- 



electron per site), with the resulting spin-orbital ex- 
change Hamiltonian for each pair {ij) : 



%ij = PijCTi ■ (Tj + Qija^aj + Rij. (8) 

Here cr^ is the spin operator at site i. The orbital opera- 
tors Ti appear in the complex form. 



(9) 



tor f? acquires an overall Zeeman field that favors ferro- 
orbital order: (r/fj) > 0. With the additional conditions 
(rf fJ) = (rf fJ) = one can prove the orbital correla- 
tions in the original basis: {y ■ Ti x Tj) = (rf f|) sin 26'. 
In the case when the {ti/t)"^ term dominates one ex- 
pects ferro-spin order (cr^ ■ aj) + \ > with easy-plane 
anisotropy for the orbital r^: (fffj) = 0, (fj^f~) ^ 0, 
implying {y ■ Ti x Tj) = {\/2){f^f~) sin2^. The ground 
state of Eq. (11) is thus expected to show strong chiral 
orbital order (y-TiXTj) ^0 throughout the wide param- 
eter range of t2/t. A full quantum-mechanical evaluation 
of the phase diagram will be carried out in a future pub- 
lication. 

Unlike the spin-DM interaction in materials, the gov- 
erning factor 7 responsible for the orbital analogue of DM 
exchange can imposed externally by the electric field in 
a controlled manner. Interesting quantum-orbital phases 
and transitions between them may be observed in a thin- 
film multi-orbital system subject to perpendicular elec- 
tric field. 
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